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Under certain circumstances, three or more interacting particles may form 
bound states. While the general few-body problem is not analytically solv- 
able, the so-called Efimov trimers appear for a system of three particles with 
resonant two-body interactions. The binding energies of these trimers are pre- 
dicted to be universally connected to each other, independent of the micro- 
scopic details of the interaction. By exploiting a Feshbach resonance to widely 
tune the interactions between trapped ultracold lithium atoms, we find evi- 
dence for two universally connected Efimov trimers and their associated four- 
body bound states. A total of eleven precisely determined three- and four-body 
features are found in the inelastic loss spectrum. Their relative locations on ei- 
ther side of the resonance agree well with universal theory, while a systematic 
deviation from universality is found when comparing features across the reso- 
nance. 
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One of the most remarkable few-body phenomena is the universally-connected series of 
three-body bounds states first predicted by Efimov (1). Efimov showed that three particles can 
bind in the presence of resonant two-body interactions, even in circumstances where any two of 
the particles are unable to bind. When the two-body scattering length a is much larger than the 
range of the interaction potential r , the three-body physics becomes independent of the details 
of the short-range interaction. Surprisingly, if one three-body bound state exists, then another 
can be found by increasing a by a universal scaling factor, and so on, resulting in an infinite 
number of trimer states (2). Universality is expected to persist with the addition of a fourth 
particle (3-7), with two four-body states associated with each trimer (5, 7); intimately tied to 
the three-body state, these tetramers do not require any additional parameters to describe their 
properties. 

Ultracold atoms are ideal systems for exploring these weakly bound few-body states be- 
cause of their inherent sensitivity to low-energy phenomena, as well as the ability afforded by 
Feshbach resonances to continuously tune the interatomic interactions. Pioneering experiments 
with trapped, ultracold atoms have obtained signatures of individual Efimov states (8-12), as 
well as two successive Efimov states (13, 14), via their effect on inelastic collisions that lead to 
trap loss. Evidence of tetramer states associated with the trimers has also been found (13, 15). 
Although the locations of successive features are consistent with the predicted universal scaling, 
systematic deviations as large as 60% were observed, and attributed to non-universal short-range 
physics (13). In the work presented here, we use a Feshbach resonance in 7 Li for which a/r 
can be tuned over a range spanning 3 decades (16). This enables the observation of multiple 
features which are compared to universal theory. 

We confine 7 Li in the \F = l,m F = 1) hyperfine state in an elongated, cylindrically symmet- 
ric, hybrid magnetic plus optical dipole trap, as described previously (16). A set of Helmholtz 
coils provides an axially oriented magnetic bias field used to tune the two-body scattering length 
a via a Feshbach resonance located near 737 G (17). For a > 0, efficient evaporative cooling is 
achieved by setting the bias field to 717 G, where a ~ 200 ciq (with a the Bohr radius), and re- 
ducing the optical trap intensity. Depending on the final trap depth, we create either an ultracold 
thermal gas just above the condensation temperature T c , or a Bose-Einstein condensate (BEC) 
with > 90% condensate fraction. For investigations with a < 0, we first set the field to 762 G 

where a 200 a and proceed with optical trap evaporation, which is stopped at a temperature 

T slightly above T c . In both cases the field is then adiabatically ramped to a final value of a and 
held for a variable hold time. The fraction of atoms remaining at each time is measured via in 
situ polarization phase contrast imaging (18) for clouds where the density is high, or absorption 
imaging in the case of lower densities. 

Analyzing the time evolution of the number of atoms in the trap determines the three-body 
loss coefficient L 3 (8, 13, 19) as well as the four-body loss coefficient L 4 (15). Recombination 
into a dimer is a three-body process since a third atom is needed to conserve both momentum 
and energy. For a > 0, the dimer can be weakly -bound with binding energy e = h 2 /(ma 2 ), 
where m is the atomic mass, while for a < there are only deeply-bound molecular dimers. 
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The recombination energy released in the collision is sufficient to eject all three atoms from the 
trap for a < 0, and for a > when e > U, where U is the trap depth. In the case of the BEC 
data, this latter condition holds for a < 5000 ao. Nonetheless, we assume that all three atoms are 
lost for any recombination event, because even for a larger than 5000 ao we observe rapid three- 
body loss. We ascribe this observation to a high probability for dimers to undergo vibrational 
relaxation collisions which result in kinetic energies much greater than U. Four-body processes 
proceed in a similar fashion (6, 15). 

The equation describing the dynamics of three- and four-body loss is 

- — = -^—L 3 (n 2 ) - ^—L 4 (n 3 ), (1) 
N dt 3! 4! 4X ' w 

where the brackets denote averages over the density distribution (17). For a thermal gas the 
spatial correlation coefficients g (3) and g (4) are respectively 3! and 4!, while for a BEC we set 
both to 1 (20, 21). We have verified that heating from recombination is small for our short 
observation times and therefore omit this effect in our analysis (15, 19). By fitting the time 
evolution of the number of atoms to the solution of Eq. 1 we extract L 3 and L 4 as a function 
of a. Figure SI shows the loss of atoms as a function of time in regimes where either L 3 or L 4 
dominates (17). Four-body loss is readily distinguished from three-body loss by the shape of 
the loss curve. 

Figure 1 shows the extracted values of L 3 across the Feshbach resonance, exhibiting the 
expected a 4 scaling (22, 23), but with several dips and peaks punctuating this trend. Two promi- 
nent peaks dominate the landscape for a < 0, which are labeled a[ and in Fig. 1A. We 
attribute these peaks to the crossings of the energies of the first two trimer states with the free 
atom threshold, thus providing additional pathways into deeply-bound molecular states (23). 
For a > 0, the dominant features are dips, indicated in Fig. 1A as a\ and a\ , corresponding to 
recombination minima. These minima are associated with the merging of the same two trimer 
states into the atom-dimer continuum, and have been attributed to destructive interference be- 
tween two different decay pathways into weakly-bound dimers (22, 23). We fit the data to 
L 3 (a) = 3C(a)ha 4 /m, where C(a) is a logarithmically periodic function characterizing effects 
from the Efimov states (17). The analytic expression for C(a) contains the location of one uni- 
versal trimer resonance a<0or recombination minimum a + > 0, and an inelasticity parameter 
r] related to the lifetime of the Efimov state (2). The observed features are fit individually to 
extract these parameters (Table 1). The universal theory describing Efimov physics (2) predicts 
a logarithmic spacing in the two-body scattering length between trimer states of e nls ° « 22.7, 
where s = 1.00624 is a universal parameter (i). Table 2 shows that the ratios a\ja\ and a~/a~ 
agree well with the universal theory. 

A local maximum in L 3 , indicated as a* 2 and shown in detail in Fig. IB, can be discerned 
between the two recombination minima a\ and a\. We associate this feature with an atom- 
dimer resonance, given its location with respect to the nearby minima. A simple model (13) 
has been proposed to explain the enhanced losses present at the atom-dimer resonance. This 
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model describes an avalanche process whereby a single dimer travelling through a collisionally 
thick gas shares its kinetic energy with multiple atoms, thereby increasing from 3 the effective 
number of atoms lost for each dimer formed (24). 

For a < 0, L 3 achieves its maximum value of ~10~ 19 cm 6 /s at a 2 . This value is reasonably 
consistent with the expected unitarity limit (19, 25). At even larger values of \a\, L 3 saturates to 
a value below the unitarity limit, a behavior previously seen in experiment (8) and in numerical 
calculations (25, 26). 

The four-body loss coefficient L 4 for a < was also extracted from the data, and the results 
are presented in Fig. 2. Three resonant peaks in L 4 are observed, which we associate with 
the crossings of tetramer states with the free atom continuum (3-7, 13, 15, 27). Two universal 
tetramers are predicted to accompany each Efimov trimer (5, 7). The solid line in Fig. 2 is 
calculated using only the observed three-body locations and widths in addition to an overall 
scaling, without any other free parameters (17). The agreement between this curve and the 
data lead us to assign the peaks to the second tetramer of the first Efimov trimer a[ 2 and both 
tetramers of the second Efimov trimer a 2l and a T 22 (15). While we do not have the resolution to 
detect an enhancement in L 4 at the expected location of the first tetramer a T xv an enhancement of 
L 3 is observed at the expected location (Fig. 1A) which we tentatively identify with a T yi (7, 13). 
The existence of two tetramer states tied to a single trimer state has also been verified in 133 Cs 
(15) and 39 K (13). 

Two additional peaks in L 3 are observed on the a > side of the resonance (Figs. 1C 
and ID). Features at these relative positions have not been previously observed or predicted, 
although they occur very close to where the two tetramer states associated with the second 
trimer are expected to merge with the dimer-dimer continuum (28). We have no explanation of 
how a dimer-dimer resonance would affect the inelastic loss rate, as we expect the dimer fraction 
to be small and consequently, the probability of dimer-dimer collisions to be negligible. One 
possibility is that they arise because of an interference effect, similar to that occurring in the 
three-body process at a\ and a 2 . Presently, we tentatively associate these features with dimer- 
dimer resonances located at a* 2 i and a* 22 . 

In Table 2 we present the relative spacings of observed loss features along with those pre- 
dicted by the universal theory. Universal scaling is expected when \a\ » r , where r is the 
van der Waals radius (33 a for Li) (29). Another requirement for universality is that \a\ » \R e \, 
where R e is the effective range (14). Figure S4 shows that R e is relatively small over the rele- 
vant field range, and is — 10ao on resonance (17). For comparison, in the |1,0) state of 7 Li, 

R e 30 ao at the resonance near 894 G (14). Both conditions for universality are well-satisfied 

for the second Efimov state, but the requirement that \a\ » r is only marginally satisfied for the 
first. Nonetheless, we find good agreement with the universal scaling relations between features 
on each side of the Feshbach resonance separately. 

The relationships of features across a Feshbach resonance are also thought to be universally 
connected (2, 26). However, when we compare features across the Feshbach resonance, we 
find a systematic discrepancy with theory of a factor of 2 (Table 2). This discrepancy can be 
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expressed as a difference in the three-body short-range phase between the two sides of the Fes- 
hbach resonance AO = s ln(\a~\/a + ) (22, 26). The locations of the features reported here result 
in phase differences of 0.92(10)(0) and 0.86(4)(17) (the uncertainties are defined in Table 1) for 
the first and second trimer, respectively, whereas the universal prediction is 1.61(3) (2). One 
of the effects of finite temperature is to both broaden the trimer resonances and to push them 
towards smaller \a\ (8, 25, 26). This would decrease the values of AO since we extract L 3 from 
a thermal cloud at tr and a much colder BEC at a + . Measurements of 39 K also show a dis- 
crepancy with theory across the resonance, but with AO = 1.91(7) (13). On the other hand, 
measurements of the first trimer resonance and second trimer recombination minimum in the 
|1, 0) state of 7 Li result in AO = 1.7(2) in good agreement with universal theory, assuming the 
universal scaling of 22.7 between trimer states (14). These variations in AO may indicate the 
need for additional physics to be included in the universal model (26, 30). 
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a > 



a < 



a\ = 119(H)(0) 



a. = -298(10)(1) 



a\ = 2676(67)(128) a~ = -6301(264)(740) 



a* = 608(H)(7) 
f 2X * 1470(15)(38)] 
2 * 3910(60)(278)] 
t]\ = 0.079(32)(20) 
77+ = 0.039(4)(10) 



'i.i 



a l2 
i 

T 



a 



2.2 



-120(20)(0) 
-295(35)(1) 
-2950(200)(150) 
-6150(800)(700) 



17- = 0.13(1)(3) 



Table 1. Locations (in a ) of three- and 
four-body loss features and inelasticity pa- 
rameters (dimensionless) (17). The fea- 
tures a* 2 j and a* 2 2 are tentatively assigned. 
The first number in parentheses character- 
izes the range over which x 2 of the fit to 
theory increases by one, while simultane- 
ously adjusting the other parameters in the 
fit. The second number characterizes the 
systematic uncertainties in the determina- 
tion of a (17). 



a > 



a < 



a — > +00 



Ratio 


Data 


Theory 


A(%) 


a 2 /a x 


22.5(22)(11) 


22.7* 


-1(9X5) 


^2 / ^2 


4.40(14)(16) 


4.46* 


-1(3)(4) 


^2 1 ^2 


* 2.42(5)(4) 


2.37* 


+2(2)(2) 


2 / ^2 


* 6.4(2)(4) 


6.6* 


-3(2)(6) 


a 2 /a~^ 


21.1(1 1)(24) 


22.7* 


-7(5)(11) 


a\ x la\ 


~0.40(7)(0) 


0.43* 


-6(16)(0) 


& j 2 1 j 


« 0.99(12)(0) 


0.90* 


+10(14)(0) 


^2 1 ^^2 


« 0.47(4)(4) 


0.43* 


+9(9)(9) 


^?2 2 


* 0.98(13)(1) 


0.90* 


+8(14)(1) 




2.5(2)(0) 


4.9* 


-49(5)(0) 


\ a 2 \/ a 2 


2.4(1X4) 


4.9* 


-52(2)(9) 


\a-\/a* 2 


0.49(2)(1) 


0.97* 


-49(2)(1) 




10.4(5)(14) 


22.0* 


-53(2)(6) 



Table 2. Relative locations of 
loss features, those predicted by 
theory, and the percent difference 
A = (Data/Theory - 1). The un- 
certainties are those propagated 
from Table 1 . 



References: * (2); f (7); * (28) 



1 



400 600 800 1000 2000 3000 3000 4000 5000 6000 

Scattering length a (a Q ) 



Fig. 1. (A) L 3 as a function of a. Data shown with (♦) correspond to a thermal gas with N ~ 10 6 , 
T ~ 1-3 /iK (31), and U ~ 6/uK and were taken with radial and axial trapping frequencies 
a> r = (2n) 820 Hz and u> z = (2^) 7.3 Hz, respectively. The remaining data correspond to a 
BEC with N ~ 4 x 10 5 , T < 0.5 T c , U ~ 0.5//K, and oj r = (In) 236 Hz. We adjust a> z (17) 
to enhance or reduce three-body loss, where co z = (In) 1.6 Hz (a), co z = (27r)4.6Hz (•), and 
a> z = (In) 16 Hz (■). The dashed lines show an a 4 scaling. The solid thick lines are fits to an 
analytic theory (2, 17). The thin green lines show the square of the energies, in arbitrary units, 
of the first and second Efimov states as predicted from the universal theory (2) where we have 
fixed the location of the first Efimov state to overlap with a\, and the atom-dimer continuum 
is coincident with the dashed line for a > 0. Several representative error bars are shown (17). 
(B — D) Detail around the loss features associated with the atom-dimer and two possible dimer- 
dimer resonances. The dotted lines are the fit to Eq. S4, while the solid lines include additional 
superimposed Gaussian fits to account for the features not described by Eq. S4. 
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Scattering length a (a ) 

Fig. 2. L4 extracted from a thermal gas. The solid curve is motivated by theory (17, 27), and the 
dashed curve is the solid curve divided by a 1 (6). The uncertainty in L 4 from the fit is a factor 
of 2, while the sytematic uncertainty is a factor of 3 due to uncertainties in co r , to z , N and T . For 
\a\ > 2 x 10 4 ao differentiation between three- and four-body losses becomes unreliable due to 
the very fast decay rates. Data with L 4 < 10" 36 cm 9 /s are consistent with no four-body loss. 
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Materials and Methods 

A set of non-Helmholtz coils are used to add or subtract additional axial confinement in the 
hybrid magnetic plus optical dipole trap used in the experiment. The radial trapping frequency 
u> r is determined from atom loss by parametric excitation, and the axial trapping frequency oj z 
is determined from collective dipole oscillations. 

Determination of Scattering Length 

The 5-wave scattering length a is controlled via a magnetic Feshbach resonance (SI). We extract 
a (for a > 0) as a function of magnetic field B from the axial size of a Bose-Einstein condensate 
(52). The measured functional form of a vs. B is well described by a Feshbach resonance 
fit a(B) = a BG [\ + A/(B - B M )], where the values a BG = -24.5^;° a , A = 192.3(3) G, and 
Boo = 736.8(2) G were previously reported (52). The standard deviation of the residuals from 
the Feshbach resonance fit is 15% for a < 10 3 ao and 30% for a > 10 3 ao (Fig. S2). 

To repeatably achieve very large values of a it is necessary to have both high field stability 
and accurate knowledge of the location of fi M . We determine the shot-to- shot stability and cali- 
bration of the magnetic field from radio frequency spectroscopy on the |1, 1) — > |2, 2) transition. 
We have improved the control of the current in the coils that provide the magnetic bias field in 
our experiment such that a Lorentzian characterizing the shot-to- shot field stability has a full 
width at half maximum of 1 15 kHz, corresponding to 42 mG at a bias field of 717 G (Fig. S3C). 
With this improved field stability we have increased the precision in the determination of the 
resonance location to B^ = 736.97(7) G. The uncertainty in is dominated by systematic 
uncertainty in the extracted values of a from the measured axial sizes (52). The fractional un- 
certainty in the determination of a is given by Sa/a = 5B/(B - B^) « 1 .5 x 10~ 5 a/ao, where 6B 
is dominated by the uncertainty in B^. 

Since we have only measured a for a > 0, we have no direct knowledge of a < 0. However, 
a coupled-channels calculation (S3) agrees with the Feshbach resonance fit to within 10% over 
the range of 10 < a/ao < 4 x 10 4 (Fig. S3) which gives us confidence that the Feshbach 
resonance fit is equally reliable on the a < side of the resonance. 
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Determination of the Loss Coefficients 



Extraction of L 3 and L 4 from the measured atom number loss curves N(t) requires the evaluation 
of the spatially-averaged moments of the density distribution (n 2 ) and (n 3 ). By comparing the 
measured distributions with a Thomas-Fermi inverted parabola in the case of a pure Bose- 
Einstein condensate, we find to a good approximation that the distributions remain in thermal 
equilibrium throughout the decay process. For a condensate, the axial Thomas-Fermi radius is 
R = (l5Ti 2 co 2 Na/m 2 co 4 z ) l/5 , the peak density is n = (15No) 2 r )/(SnR 3 o) 2 z ), and (n 2 ) = y 2 ' 5 N 4/5 , 
where y = (25 m 6 u) 4 u) 2 )/ (627 2 y/42n 5 h 6 a 3 ). The observed decay fits well to a purely three-body 
loss process for a condensate, so we neglect L 4 in this case. Since we are not explicitly fitting 
for L 4 , four-body effects if present may lead to an increase in the extracted loss rate L 3 (S4). 
The decay is then described by 



which has the solution 



^ :s74- (S3) 



A thermal gas is well described by a cylindrically- symmetric Gaussian where (n 2 ) = n 2 / V27, 
(n 3 ) = n 3 /S, and the peak density is n p = N(a> z / aj^imoj 2 . /2nkBT] 3/2 . Heating due to recom- 
bination is expected to become important when e < U (S5). However, there is no appreciable 
change observed in the Gaussian width during the decay even though the loss mechanism pref- 
erentially targets atoms at higher densities. This may be due to a lack of rethermalization during 
the decay (S6). We find that both L 3 and L 4 contribute to the loss for the thermal gas. Since we 
have not found a closed-form solution to Eq. 1, we instead use the following implicit solution 
to extract L 3 and L 4 : 



3V3 
2^1," 



(N \ 2 1 27L 4 / N \ 81V3/2 \/ N \ 

\n) - 1 \ + s^ 2 V-^)--6jir log Wo) 



N \ 8 V3L 3 + 9L 4 n p 



kNq ) 8 V3L 3 + 9L 4 nJN /N) 



(S4) 



where we have assumed g (3) = 3! and g (4) = 4! for a non-condensed gas. 

In Fig. 1 the vertical error bars correspond to the range in L 3 for which the;^ 2 of the fit to 
Eq. S3 increases by one, while simultaneously adjusting L 4 and Af to minimize^ 2 . Systematic 
uncertainties in oj r , co z , N, and T, which are not included in these error bars, contribute as much 
as a factor of 2 in the uncertainty of L 3 . The representative horizontal error bars are due to shot- 
to- shot variation in the magnetic field and the determination of a from the Feshbach resonance 
fit. Background loss limits the sensitivity of the measurement to L 3 > 2(1) x 10" 28 cm 6 /s. The 
error bars in Fig. 2 are similarly determined. 
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Comparing with Theory 

The universal theory (57) describing Efimov physics predicts that the three-body loss rate coeffi- 
cient is described by L 3 (a) = 3C{a)Tia A I m where C(a) is a logarithmically periodic modulation. 
The following expression describes this modulation: 



C(a) = 



4590 sinh(27/-) 



sin (so\n(a/a )) + sinh rj 
61.12e- 2 " + [sin 2 (j ln(a/a + )) + sinhV] + 16.84(1 - e~ 4ri+ ) {a > 0) 



(a < 0), 

(S5) 



where the first and second terms for a > account for coupling to weakly- and deeply-bound 
dimer states, respectively (57, S8). The value <sT denotes the resonance location when the energy 
of the Efimov trimer is degenerate with the free atom continuum, and the value a + is the location 
of a recombination minimum (59). This expression is log-periodic with C(e n/s °a) = C(a), where 
the universal parameter so = 1.00624 is known from theory (57, S10). 

The four-body loss coefficient L\ is predicted to have a similar form to that of L 3 : 

j , t, An ^M 7 sinh(2?T) 

L 4 (a,a ) = 4C 4 = = — (a < 0), (S6) 

m sin (s ln(a/a T )) + sinh rf 

where C4 is a theoretically undetermined universal constant (Sll). Eq. S5 is phenomenologi- 
cally derived from the theory of Ref. Sll (572). We find that C 4 = 16(8) x 10 4 in the region 
1000 < -a/ao < 2500, assuming that 77" = 0.13, as for the three-body resonance. In Fig. 2 we 
plot \ {L 4 (a, 0.90 a[) + L 4 (a, 0.43 ap} where we have replaced a T with the predicted locations 
of the two tetramer states linked to the first trimer state (S4, Sll). 
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Fig. SI. Loss dynamics at two values of a < for a thermal gas. The dots are data. The dotted 
red line is a fit of the data to the solution of Eq. 1 with only three-body loss accounted for, the 
dashed blue line is the fit when only four-body loss is included, and the solid green line is a 
fit accounting for both effects (Eq. S3). (A) a = -1800 ao, where three-body losses dominate; 
(B) a = -3300 ao, near a\ x where four-body losses dominate. 
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Fig. S2. (A) a extracted from the axial size of Bose-Einstein condensates as a function of 
magnetic field. Results of a coupled-channels calculation are shown by the solid red line. The 
dashed black line is the Feshbach resonance fit. (♦) Data previously reported with trapping 
frequencies to r = (2n) 193 Hz and a> z = (2n) 3 Hz (52). Data with a> r = (2n) 236 Hz and 
a> z = (2n) 4.6 Hz (•) or co z = (In) 16 Hz (■). Beyond mean field effects become important when 
n a 3 > 0.1 (SI 3). We apply a mean field correction for data with 0.1 < n a 3 < 1, and omit data 
with n a 3 > 1 in the Feshbach resonance fit (52). (B) Full range of data spanning 7 decades 
in a. (C) Fractional residuals of the extracted values of a from the Feshbach resonance fit. 
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Fig. S3. (A) a vs. magnetic field from a coupled-channels calculation. (B) Fractional difference 
between the coupled-channels calculation and the Feshbach resonance fit used to determine a 
(solid red line a > 0, dashed blue line a < 0). (C) Radio frequency spectroscopy signal at 717 G 
showing a full width at half maximum of 1 15 kHz. 
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Fig. S4. The effective range R e (solid red) and scattering length a (dashed blue) vs. magnetic 
field, extracted from a coupled-channels calculation through a low energy expansion k cot 6 = 
-l/a + R e k 2 /2, where 6 is the scattering phase shift (SI). The dotted vertical line is the location 
of B m . 
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